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ABSTRACT 

Currently popular models for progenitors of gamma-ray bursts (GRBs) are 
the mergers of compact objects and the explosions of massive stars. These two 
cases have distinctive environments for GRBs: compact object mergers occur in 
the interstellar medium (ISM) and the explosions of massive stars occur in the 
preburst stellar wind. We here discuss neutrino afterglows from reverse shocks 
as a result of the interaction of relativistic fireballs with their surrounding 
wind matter. After comparing with the analytical result of Waxman & Bahcall 
(2000) for the homogeneous ISM case, we find that the differential spectrum 
of neutrinos with energy from ~ 3 x 10 15 to ~ 3 x 10 17 eV in the wind case 
is softer by one power of the energy than in the ISM case. Furthermore, the 
expected flux of upward moving muons produced by neutrino interactions below 
a detector on the surface of the Earth in the wind case is ~ 5 events per year 
per km 2 , which is about one order of magnitude larger than in the ISM case. In 
addition, these properties are independent of whether the fireballs are isotropic 
or beamed. Therefore, neutrino afterglows, if detected, may provide a way of 
distinguishing between GRB progenitor models based on the differential spectra 
of neutrinos and their event rates in a detector. 



Subject headings: gamma-rays: bursts - elementary particles - shock waves - 
stars: mass loss 



1. Introduction 

The study of gamma-ray bursts (GRBs) has been revolutionized due to observations 
of multiwavelength afterglows in the past few years, but the nature of their progenitors 
remains unknown (for a review see Meszaros 1999). Two currently popular models for 
GRB progenitors are the mergers of compact objects (neutron stars or black holes) and 
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the explosions of massive stars. In the former model, compact objects are expected 
to have significant spacial velocities so that their mergers would take place at many 
kiloparsecs outside their birthplaces. Thus, GRBs produced by this model would occur in 
the interstellar medium (ISM) with density n ~ 1 cm -3 . Strong evidence for the massive 
star progenitor model has been recently discovered. GRB 980425 was probably associated 
with the relatively nearby Type Ic supernova (SN) 1998bw (Iwamoto et al. 1998; Kulkarni 
et al. 1998), and the supernova-like emission was also found in GRB 980326 (Bloom et al. 
1999) and GRB 970228 (Reichart 1999; Galama et al. 2000). These observations show that 
some or possibly all long-duration GRBs arise from the core collapse of massive stars. It 
has been widely believed that GRBs associated with supernovae should unavoidably occur 
in the preburst stellar wind environment with mass density p oc R~ 2 . If GRB emission is 
isotropic, the X-ray and optical afterglow in the wind case must decline more rapidly than 
in the ISM studied analytically by Dai & Lu (1998), Meszaros, Rees & Wijers 

(1998), Panaitescu, Meszaros & Rees (1998) and Chevalier & Li (1999, 2000). Guided by 
this argument, Dai & Lu (1998) suggested for the first time that GRB 970616 is a possible 
wind interactor based on the rapid fading indicated by two X-ray flux measurements. 
Recently Chevalier & Li (1999) argued that GRB 980519 is an excellent wind interactor 
based on its X-ray, optical and radio data (although these observational data were shown 
analytically and numerically to be consistent with the dense medium model by Dai & Lu 
[2000] and Wang, Dai & Lu [2000]). Furthermore, the afterglow data of some other bursts 
(e.g., GRB 970228, GRB 970508, GRB 980326 and GRB 980425) are also consistent with 
the wind environment model (Chevalier & Li 2000). These afterglows have been argued to 
be further evidence for massive stars as GRB progenitors. 

In this paper we study neutrino afterglows from reverse shocks as a result of the 
interaction of relativistic fireballs with their surrounding wind matter by assuming that 
GRBs result from the explosions of massive stars. We find that the differential spectrum 
of neutrinos below ~ 3 x 10 15 eV is proportional to e^ 1 but the differential spectrum of 
neutrinos with energy from ~ 3 x 10 15 to ~ 3 x 10 17 eV steepens by one power of the 
energy. In addition, the expected flux of upward moving muons produced by neutrino 
interactions below a detector on the surface of the Earth is ~ 5 events per year per km 2 
for typical parameters. We also find that this flux is ~ 10 times larger than estimated by 
Waxman & Bahcall (2000), who studied neutrino emission from reverse shocks produced by 
the interaction of fireballs with the interstellar medium (ISM) . 

The neutrinos are produced by 7r + created in interactions between accelerated protons 
and synchrotron photons from accelerated electrons in a relativistic fireball. This neutrino 
emission during the GRB phase was studied in the internal shock models by Waxman & 
Bahcall (1997), Halzen (1998) and Rachen & Meszaros (1998). It was found that a fraction, 
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~ 0.1, of the fireball energy would be converted by photomeson production to a burst of 
neutrinos with typical energy of a few 10 14 eV (but also see Vietri 1998). The property 
of such neutrino bursts is independent of whether the ambient matter is a stellar wind 
or a constant density medium. Observations of these bursts could be used to test the 
simultaneity of neutrino and photon arrival to an accuracy of ~ 1 s, the weak equivalence 
principle, and the vacuum neutrino oscillation theory (Waxman & Bahcall 1997). 

The structure of this paper is as follows: In Section 2 we analyze reverse shocks 
produced during the interaction of ultra-relativistic fireballs with the surrounding wind 
matter and discuss the photon emission from these shocks. In Section 3 we investigate 
neutrino afterglow emission as a result of photo-meson interaction in the reverse shocks 
and in Section 4 we discuss the detectability of such afterglows. In the final section, several 
conclusions are drawn. 



2. Shock Model and Photon Emission 

We first assume that a relativistic GRB shell will interact with the surrounding stellar 
wind via two shocks: a reverse shock and a forward shock. The forward shock runs forward 
into the wind while the reverse shock sweeps up the shell material. The recently observed 
optical flash of GRB 990123 has been argued to come from a reverse shock (Akerlof et 
al. 1999; Sari & Piran 1999; Meszaros & Rees 1999). We believe that the reverse shock 
emission should be common for all GRBs. The shocked ambient and shell materials are 
in pressure balance and are separated by a contact discontinuity. We assume that these 
shocked materials are uniform and move together. Sari & Piran (1995) and Mitra (1998) 
considered the jump conditions for relativistic shocks and found the common Lorentz factor 
of the shocked materials measured in the unshocked medium frame, 

£l/4pl/2 

7 = 



y/2 ' 

where Y is the Lorentz factor of the unshocked shell measured in this frame and £ = p sh / p w 
is the ratio of proper mass densities of the unshocked shell and the unshocked ambient 
medium. The proper mass density of the ambient medium is expressed as 

where M w and V w are the mass loss rate and wind velocity of the progenitor star, and 
A = M w /(AttV w ) = lO- 5 M yr-7(47r x 10 3 kms- 1 )A* = 5 x 10 n gcm-M* and R is the 
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radius of the shell in units of 1 cm (Chevalier & Li 1999, 2000). The proper mass density of 
the unshocked shell is given by 

- E ° (i\ 
Psh ~ 47ri? 2 r 2 c 2 A' [6) 

where Eq and A are the energy and the width (measured in the unshocked medium frame) 

of the initial shell. A typical value of A* ~ 1 for Wolf-Rayet stars is found from stellar 

mass- loss rates and wind velocities (Willis 1991; Chevalier & Li 2000 ). Since GRBs are 

believed to come from internal shocks, A is approximately equal to the speed of light times 

the GRB durations and thus its typical value should be ~ 10 light seconds. The rapid 

variability of GRBs and their nonthermal spectra require that T be a few hundreds (Woods 

& Loeb 1995). From the observed fluences of some GRBs and their measured redshifts, Eq 

is estimated to be between 10 52 and 10 54 ergs. A recent analysis by Freedman & Waxman 

(2000) also gives this estimate. Scaling the involved quantities with these typical values, we 

find 

^ = 655-^^, (4) 

where -E53 = i?o/10 53 ergs, r 30 o = T/300 and Ai is in units of 10 light seconds. From 
equations (1) and (4), the Lorentz factor of the shocked shell material is rewritten as 



E, 



1/4 

53 



7 = 62 <5) 

Following Sari & Piran (1995) and Mitra (1998), we further derive the Lorentz factor of the 
shocked shell material measured in the unshocked shell rest frame, 

|^ = ir = 242 ^r 3W (6) 

V2 2 7 EH 1 

which implies that the reverse shock is relativistic. After the reverse shock passes through 
the shell, the shock front disappears. Instead of maintaining a constant Lorentz factor (e.g., 
equation [6]), the shocked materials slow down with time based on the Blandford-McKee 
(1976) self-similar solution. In the following we discuss photon emission from the reverse 
shock. 

Because of pressure balance across the contact discontinuity, the shocked shell material 
and the shocked wind material have not only the same bulk Lotentz factor but also the 
same internal energy density. According to relativistic shock jump conditions, we obtain 
the internal energy density of the shocked shell material, 



e' = 4 7 2 p w c 2 = 2r Pw c 2 e 1/2 . 



(7) 
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We assume that e e and e B are the fractions of the internal energy density (in the shocked 
material rest frame) that are carried by electrons and magnetic fields respectively. 
The minimum Lorentz factor of the electrons accelerated behind the reverse shock is 
approximated by 7 m ~ {m, p /m e )e e ^' (Waxman & Bahcall 2000), viz., 

' 4 l/4 A l/4 r 

A A r I A * ^10 1 300 | / c x 

7 m « 445e e ,_i —j-. , (8) 



E, 



1/4 
53 



where e e _i = e e /0.1. Moreover, the magnetic field strength in the shocked shell material is 
given by 

h 

where U = R/c is the time in the burster's rest frame. Substituting the relation between 
this time and the observed time (t), tj, = 2^ 2 t/{l + z), into the above equation, and using 
equations (4) and (5), we further have 

= 3.8 x 10^.3 (i±£)(^!) G , (10) 

\ 53 5 J 

where e B -s = e B /10~ 3 , z is the redshift of the source and t s — t/1 s. It should be noted that 
e e ~ 0.1 (Freedman & Waxman 2000), but €b is highly uncertain and its reasonable value 
may be taken from ~ 10~ 2 to ~ 10 -6 . Several previous studies of GRB afterglows (e.g., 
Galama et al. 1999; Dai & Lu 1999, 2000; Wang, Dai & Lu 2000) gives e B ~ 10~ 4 - 10" 6 . A 
recent detailed study of the afterglows of GRBs 980703, 990123 and 990510 by Panaitescu 
k Kumar (2000) leads to e B ~ 10~ 3 - 10~ 4 . In addition, Holland et al. (2000) find that e B 
is as small as 10~ 5 . Therefore, we choose a typical value: e B ~ 10~ 3 . 

Let's consider synchrotron radiation of the electrons accelerated behind the reverse 
shock. We first derive two characteristic frequencies of synchrotron photons: the typical 
frequency v m corresponding to the minimum electron Lorentz factor and the cooling 
frequency v c . From equations (5), (8) and (10), we obtain the typical frequency in the 
observer frame, 



m = — — ^ — — — — = 6.3 x 10 16 £e 2 ,ef 3 ^gUgo Hz (11) 
l + z2nm e c e,-i b,-s 



The cooling frequency corresponds to the cooling Lorentz factor j c , at which an electron 
is cooling on the dynamical expansion time. We believe that this Lorentz factor in the 
reverse shock will increase with time because of the cooling timescale oc B'~ 2 oc t 2 . 
Initially, 7 C — 1 < 1. At this stage, a cooling electron may be non-relativistic and its 
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kinetic energy E e « (l/2)^m e c 2 f3 2 . Thus, its cooling timescale due to cyclotron radiation, 
measured in the observer's frame, can be estimated as to — (1 + z)E e / P cyc ((3), where 
Pcyc(P) = (4/3)o"tC7 2 /9 2 -B /2 /(87t) is the cyclotron radiation power (in the local observer 
frame) of an electron with velocity of f3c and with <jt being the Thomson scattering cross 
section. Using equation (10), we easily find 

^'■^tP^ 1 < 12 > 

-^53 

Please note that to is independent of (5. This implies that at t < t an electron accelerated to 
7 m in the magnetic field B' will be able to cool to become non-relativistic, initially through 
synchrotron radiation and subseqently through cyclotron radiation, on the dynamical 
expansion time t. However, when t > to, the magnetic field {B 1 oc t~ v ) will become weaker 
and the cooling timescale due to cyclotron radiation must be longer than t, so an electron 
with 7 m cannot cool to a non-relativistic velocity on time t. In this case, cyclotron radiation 
is no longer a cooling mechanism but it should be replaced by synchrotron radiation. Since 
for typical parameters t ~ 1 s, which is much less than the durations of long GRBs from 
the collapse of massive stars, we will discuss the photon and neutrino emission from the 
reversely shocked matter at t > to in the remaining text. According to Sari, Piran & 
Narayan (1998), the cooling Lorentz factor is defined by 

77 c m e c 2 = P syn {j c )t/ (1 + z), (13) 

where P sy n(7c) = (4/3)<ttC7 2 7 2 /3 2 -B /2 /(87t) is the synchrotron radiation power (in the local 
observer frame) of an electron with Lorentz factor of 7 C (and f3 ~ 1). Equation (13) leads to 

67rm e c(l + z) , f l + z y 1 E^t a 

7c= *ryB»t =2 ^-A—) a^W- (14) 

It is clear that 7 C 3> 1 for t 3> 1 s, implying the cooling electrons are indeed relativistic. 
Using this equation, we further derive the cooling frequency in the observer frame, 

- = ,3 x (I±£y ^ Hz. (15) 



1 + z 27im e c ' B >~ s V 2 J Al A^ 2 



We can see from equations (11) and (15) that for typical parameters the cooling frequency 
is much lower than the typical frequency, indicating that all radiating electrons cool rapidly 
down to the cooling Lorentz factor. In other words, the shocked shell material is in the 
fast cooling regime. It is interesting to note that this conclusion has also been drawn by 
Chevalier & Li (2000). Therefore, the observed specific luminosity peaks at e c = hu c rather 
than e m = hu m , with a peak value approximated by 

1/2 

L ec = (2nh)-\l + z) lNe P' c = 2.7 x 10 62 e^ 3 (^) s~\ (16) 

V 1 / 1 300^10 
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where N e = E ct/[(1 + z)Tm p c 2 A] is the number of radiating electrons in the shocked shell 
region, and P' c = m e c 2 a T B' /(3e) is the power radiated per electron per unit frequency in 
the shocked shell rest frame. 

We turn to derive the synchrotron self-absorption frequency of the reversely shocked 
matter. In the comoving frame of the shocked matter, the absorption coefficient for u' > v' c 
is given by 

a > = ^! ( —^—Y /2 (m^r^KXB'^y'u'-^^T (*P±*) r f + 
" 8nm e \2nm 3 e c 5 J K J V 12 / V 12 / ' v ; 

where K = 4(p— l)7'7 c (p s h/m p ) and A = (1/2) J 7r (sina)( p+2 )/ 2 smada (Rybicki & Lightman 
1979). In the present case, the electron distribution index p = 2 because v c u m , and the 
width of the reverse shock R' = / ~fct/(l + z). Setting T(u' a ) = a' u , R' = 1, we can derive the 
synchrotron self-absorption frequency in the obsever's frame, 

„ a ^ =2 . 2xl0 . 5 (i±£)- i/3 5g: 3 H, as) 

l + Z V 2 / A l/6 r l/3,2/3 V ! 

Hence, u c <^ u a <^ u m for typical parameters. 

We assume that the electrons behind the reverse shock follow a power law energy 
distribution, dn' e /d^ e oc % 2 for 7 e > 7 m (Blandford & Eichler 1987). In this case, the 
synchrotron radiation spectrum is a broken power law: 

^ c (e a Ac)- 1/2 (e 7 Aa) 5/2 if e 7 <e a 

£e c (e 7 /e c r 1/2 if e a <e 7 <e m (19) 

Le c {^rn/ ^c) ^ (^7/^m) if £7 ^ frai 

where e a = hu a . The protons behind the reverse shock are expected to be accelerated to the 
same power-law distribution as the electrons (with the maximum proton energy which will 
be estimated in the next section). 



3. Neutrino Emission 

For convenience, we hereafter denote the particle energy measured in the shocked shell 
rest frame with a prime, and the particle energy in the observer frame without prime, 
e.g., e 7 = 7e 7 /(l + z). We now consider neutrino production in the wind case. Assuming 
n 7 (e 7 )<ie 7 to be the photon number density in the shocked shell rest frame and following 
Waxman & Bahcall (1997), we can write the fractional energy loss rate of a proton with 
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energy e' p due to pion production, 

1 det 



i 1 dt \ r°° r°° 

K\0 = = 7r-5 c / tfa7 w (e)f(e)e / dxx- 2 n'{x), 

C p dt' Je Je/2j p 



(20) 



where 7 p = e' p /m p c 2 , a n (e) is the cross section for pion production for a photon with energy 
e in the proton rest frame, £(e) is the average fraction of energy lost to the pion, e = 0.15 
GeV is the threshold energy, and the photon number density is related to the observed 
specific luminosity by n'(x) = L eri ( / -fx)/[4irR 2 c(l + 2)72;]. Because of the A resonance, we 
find that photo-meson production is dominated by the interaction with photons in the 
energy range e 7 > e m . Considering the photon spectrum in this energy range, equation (20) 
leads to 

.'-!( r\ L £c ( e 'c \ ( lp e m \ Cpeak^peakAe 

* [p) 37r(l + z)R^\e' m ) \e peak ) e pcak ' 1 ) 

where a pea k = 5 x 10~ 28 cm" 2 and ^ pcak = 0.2 at the resonance e = e peak = 0.3 GeV, and 
Ae = 0.2 GeV is the peak width. The fraction of energy loss of a proton with observed 
energy e p by pion production, f n (e p ), is defined by t'~ l times the expansion time of the 
shocked shell material (~ R/jc). Thus, we have 

It is interesting to note that fn{^p) is independent of and T. Similarly to the cooling 
electron Lorentz factor defined by Sari et al. (1998) (see equation [13]), we can define 
the cooling proton energy e P)C based on f n (e PtC ) = 1. According to equation (22), we find 
e p c ps 5 x 10 16 eV for typical parameters. This implies that the protons with energy > e PjC 
accelerated behind the reverse shock must lose almost all of their energy (viz., significant 
cooling) due to photo-meson interactions, but the protons with energy < e PiC lose only a 
fraction (~ f n ) of their energy. 

We now turn to discuss the neutrino spectrum. The photo-meson interactions include 
(1) production of n mesons: jry — > p + ir° and pj — > n + n + , and (2) decay of n mesons: 
7T° — > 27 and 7r + — > /i + + — > e + + u e + + v^. These processes produce neutrinos with 
energy ~ 0.05e p (Waxman & Bahcall 1997). Since the protons with energy < e PiC lose only a 
fraction (~ f w oc e p ) of their energy, the differential spectrum of neutrinos below the break 
energy ~ 3 x 10 15 eV is harder than the proton spectrum by one power of the energy. But 
since the protons with energy > e PjC accelerated behind the reverse shock must lose almost 
all of their energy, the neutrino spectrum above the break traces the proton spectrum. 
Therefore, if the differential spectrum of accelerated protons is assumed to be a power law 
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form n(e p ) oc e p 2 , the differential neutrino spectrum is n{e v ) oc e v x below the break and 
n(e u ) oc e~ 2 above the break. 

The maximum energy of the resultant neutrinos is estimated as follows. This energy is 
determined by the maximum energy of the protons accelerated by the reverse shock. The 
typical Fermi acceleration time is t' a = fR L /c, where R L — (1 + z)e p /( / yeB') is the Larmor 
radius and / is of order unity (Hillas 1984). The requirement that this acceleration time 
is equal to the time for energy loss of protons (t' n ) due to pion production leads to the 
maximum proton energy, 

W = 5.6 x 10^f-y\^ef_ 3 (±±£) eV. (23) 

From this equation, we can draw two conclusions: (1) For reasonable parameters, the 
maximum proton energy is ~ 6 x 10 18 eV, which is two orders of magnitude smaller than 
the maximum energy of the protons accelerated by the reverse shock in the ISM case 
(Waxman & Bahcall 2000). The physical conditions in the reverse shock for the ISM case 
imply that protons can be Fermi accelerated to ~ 10 21 eV (Waxman 1995; Vietri 1995; 
Milgrom & Usov 1995; see Waxman [1999, 2000] for recent reviews). (2) The maximum 
energy of neutrinos produced in the wind case is ~ 3 x 10 17 eV. 



4. Detectability 

We discuss the detectability of the neutrino afterglow emission in the wind case. 
Since the protons with energy > 5 x 10 16 eV must lose almost all of their energy due 
to photo-meson interactions, the present day neutrino energy density due to GRBs is 
approximately given by U v = (l/2)(l/2)t#i?, where the first factor 1/2 accounts for the 
fact that about one half of the proton energy is lost to neutral pions which do not produce 
neutrinos, the second factor 1/2 accounts for the fact that about one half of the energy 
in charged pions is transferred to + and tu = 10 Gyr is the Hubble time. Here we 
assume that E = 0.8 x 10 44 erg Mpc~ 3 yr -1 is the production rate of GRB energy per unit 
volume (Waxman & Bahcall 2000). The neutrino flux is thus approximated by 

^^xlO-^^-'cnrW. (24) 

The resulting high-energy neutrinos may be observed by detecting the Cherenkov light 
emitted by upward moving muons produced by neutrino interactions below a detector on the 
surface of the Earth (Gaisser, Halzen & Stanev 1995; Gandhi et al. 1998). Planned 1 km 3 
detectors of high energy neutrinos include ICECUBE, ANTARES, NESTOR (Halzen 1999) 
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and NuBE (Roy, Crawford & Trattner 1999). The probability that a neutrino could produce 
a high-energy muon in the detector is approximated by P u ^^ ~ 6 x 10 _4 (e Jy /3 x 10 15 eV) ' 5 . 
Using equation (24), we obtain the observed neutrino event rate in a detector, 

N evCQts = 2^<j> v P v ^ » 5 ( 3 ~ 1015^ ) ^ km^yr" 1 . (25) 

This equation shows that a km 2 neutrino detector should detect each year about 5 neutrinos 
(with energy of ~ 3 x 10 15 eV) correlated with GRBs. For a GRB, its neutrino emission 
from the reverse shock in the wind case should be delayed to a few seconds after the main 
burst. Waxman & Bahcall (2000) have found f n ~ 0.1 for neutrino emission from reverse 
shocks in the ISM case (where the typical energy of neutrinos is ~ 3 x 10 eV). Using the 
same expression of P u ^n, we have re-derived their neutrino event rate in a detector and 
obtained N cvcnts ~ 0.5km _ yr _1 , which is smaller than our event rate by a factor of ~ 10. 



5. Discussion and Conclusions 

Neutrino bursts during the GRB phase were studied in the internal shock models by 
Waxman & Bahcall (1997) and Halzen (1998) who found that the neutrino event rate in a 
detector (mainly neutrinos with typical energy of a few 10 14 eV) is ~ 26 events per year per 
km 2 , which is larger than our event rate by a factor of ~ 5. Compared with the analytical 
result of Waxman & Bahcall (2000), our discussions on neutrino afterglows in the wind 
case can lead to the following conclusions: (1) The protons with energy > 5 x 10 16 eV 
must lose almost all of their energy due to photo-meson interactions and thus the neutrino 
afterglow emission in the wind case is dominated by neutrinos with energy ~3x 10 15 eV. 
(2) The maximum energy of the protons accelerated behind the reverse shock in the wind 
case is ~ 6 x 10 18 eV, so ultrahigh energy cosmic rays cannot be produced in this case. 
In addition, the maximum neutrino energy is ~ 3 x 10 17 eV. (3) The neutrino differential 
spectrum below ~ 3 x 10 15 eV is proportional to e^ 1 but the spectrum between ~ 3 x 10 15 
and ~ 3 x 10 17 eV steepens by one power of the energy. (4) The observed neutrino event 
rate in the wind case is ~ 10 times larger than the one in the ISM case. 

If GRB emission is isotropic, the optical afterglow in the wind case must decline more 
steeply than in the ISM case. This has been suggested as a plausible way of distinguishing 
between the GRB progenitor models (Chevalier & Li 1999, 2000). It is widely believed 
that GRBs may come from jets (Kulkarni et al. 1999; Castro- Tirado et al. 1999). As 
argued by Rhoads (1999) and Sari, Piran & Halpern (1999), the optical afterglow from a 
jet is likely to decay more rapidly at late times than at the early stage due to the lateral 
spreading effect. If this effect is true, however, both ISM and wind cases should show the 
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same emission feature during the lateral spreading phase, and in particular on a timescale 
of days the wind density is similar to typical ISM densities so that an interaction with the 
wind would give results that are not too different from the ISM case (Chevalier & Li 2000; 
Livio & Waxman 1999). If GRBs are beamed, thus, their optical afterglow emission could 
not be used to discriminate the massive star progenitor model from the compact binary 
progenitor model. However, the neutrino afterglow emission discussed here is independent 
of whether the fireballs are isotropic or highly collimated. Therefore, neutrino afterglows, if 
detected, may be used to distinguish between GRB progenitor models based on differential 
spectra of observed neutrinos and their event rates in a detector. 

What we want to point out is that the above conclusions are drawn by considering 
typical values of the parameters entering the fireball shock model. In fact, these parameters 
may have respective distributions. It is interesting to note that such distributions may lead 
to an event rate larger than estimated in this paper (Halzen & Hooper 1999). 
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that enabled us to improve our paper, Y. F. Huang for reading carefully this paper, and X. 
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Science Foundation of China (grants 19825109 and 19773007), and by the National Project 
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